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ABSTRACT 

In order to understand the formation mechanism of the disks around Be stars it is imper- 
ative to have a good overview of both the differences and similarities between normal B stars 
and the Be stars. In this paper we address binarity of the stars. In particular, we investigate a 
previous report that there may be a large population of sub-arcsecond companions to Be stars. 
These had hitherto not been detected due to a combination of a limited dynamic range and 
spatial resolution in previous studies. We present the first systematic, comparative imaging 
study of the binary properties of matched samples of B and Be stars observed using the same 
equipment. We obtained high angular resolution (0.07-0.1 arcsec) K band Adaptive Optics 
data of 40 B stars and 39 Be stars. The separations that can be probed range from 0.1 to 8 
arcsec (corresponding to 20-1000 AU), and magnitude differences up to 10 magnitudes can in 
principle be covered. After excluding a few visual binaries that are located in regions of high 
stellar density, we detect 1 1 binaries out of 37 Be targets (corresponding to a binary fraction of 
30 ± 8 %) and 10 binaries out of 36 B targets (29 ± 8%). Further tests demonstrate that the B 
and Be binary systems are not only similar in frequency but also remarkably similar in terms 
of binary separations, flux differences and mass ratios. The minimum physical separations 
probed in this study are of order 20 AU, which, combined with the similar binary fractions, 
indicates that any hypotheses invoking binary companions as responsible for the formation of 
a disk need the companions to be closer than 20 AU. Close companions are known to affect 
the circumstellar disks of Be stars, but as not all Be stars have been found to be close binaries, 
the data suggest that binarity can not be responsible for the Be phenomenon in all Be stars. 
Finally, the similarities of the binary parameters themselves also shed light on the Be forma- 
tion mechanism. They strongly suggest that the initial conditions that gave rise to B and Be 
star formation must, to all intents and purposes, be similar. This in turn indicates that the Be 
phenomenon is not the result of a different star formation mechanism. 

Key words: techniques: high angular resolution - stars: binaries - stars: emission line, Be - 
stars: statistics 



1 INTRODUCTION 

Be stars are rapidly rotating stars surrounded by an ionized circum- 
stellar disk, whose emission is best known for the distinct doubly 
peak ed Ha lines. About 15% of all B stars belong to the Be cate- 
gory jjaschek & Jasche3l99Ct chap. 9), making them a significant 
component of the hot, massive star population. However, the ori- 
gin of the circumstellar disks around Be stars is still a n unsolved 
mystery (see e.g. the review bv lPorter & Riviniusll2003l) . 

Binarity has sometimes been proposed to be r espon sible for 
the Be phenomenon. For example. lHarmanec et al.l ( 120020 investi- 



* Based on observations collected at the European Southern Observatory 
(ESO), Paranal, Chile under programme ID 075.C-0475(A). 



gate the effects of a detached, secondary component's gravitational 
field on the mass eject ion into the equatorial plane from the pri- 
mary. iGies et alj ( fl99^> discuss the rotational history of 4> Per, and 
find the (compact) companion has had a significant tidal effect on 
the central Be star. In addition, it should be noted that companions 
are also known t o affect the disk's properties such as its density law 
or even size (e.g. ljones e t al. 2008 on k Dra). 

Most models based on binarity have in common that close bi- 
naries with separations of order AU, a re responsible for the for- 
mation of a disk (for an overview, see iNegueruelal 120071) . When 
one considers the more general case of binaries, it is plausible that 
disks are affected during the star formation stages by the presence 
of a companion star, even in the case of a wide binary. However, to 
be able to address binarity as a crucial ingredient for the existence 
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Table 1. The target sample. Unless otherwise noted in the text, photometry and spectral types are taken from the Bright Star Catalogue iHoffleit & W arren jj 
1982), the coordinates are J2000. HR 2577 is an unresolved spectroscopic binary. 



of Be stars, we need to know whether Be stars are found in binary 
systems more often than normal B stars. If so, this would indicate 
that binarity must have something to do with the production of the 
conditions required for disk formation - either directly, where the 
tidal forces acting within a (close) binary system result in a disk 
around a mass losing object, or indirectly, where the influence of 
the binary companion allows the existing disk to survive. 

The most straightforward manner to determine whether the 
Be binary fraction is significantly different from that of "normal" 
B-type stars and thereby to assess whether existing binary mod- 
els are viable in the first place is to establish whether Be stars are 
more frequently found in binaries than normal B stars. The most re- 
ce nt comparative study on Be binarity was performed 25 years ago 
bv lAbt&Cardonal ( ll984l) . In this study, incorpora ting the spectro- 
scopic survey of Be and B stars bv lAbt & Levy! 19781, and including 
both spectroscopic and visual binaries, they found the Be binary 



fraction to be 30%, similar to that of B stars. However, this paper 
was essentially a compilation of the then, inevitably incomplete, 
literature on visual binaries and was never repeated by a homoge- 
neous dedicated investig ation. Later s t udies, such as the Speckle 
interferometric study by Maso n et al. I Jl997h concentrated on Be 
stars alone. Despite their high angular resolution, Mason et al. did 
not discover many new binary systems. This is most likely due to 
their limited dynamic range - only companions which were fainter 
by up to 3 magnitudes than the primary could be detected. 

The original motivation for the current study was the unex- 
pected result that in an investigation into the spectro-astrometric 
signatures of the Ha emission lines of a random sample of 
8 Be stars, 5 were found to be a bin ary with sepa rations at 
the sub-arcsecond to arcsecond scale teainesl |2004| . see also 
lOudmaiier et alj2 008). The parameter space probed by the spectro- 
astrometry (from 0.1 arcsec to 1-2 arcsec, up to 6 magnitudes dif- 
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Figure 1. Example of the data. Shown in this 3x3 arcsec image is the 
known B-type binary HR 1892. The previously known and catalogued com- 
panion is located at a distance of 1.6 arcsec towards the South- West. The 
close companion at 0.15 arcsec (double the spatial resolution in this image) 
in the North-West direction is a new detection. 



ference, iBaines et alf 2006) was not covered in any previous study 
into the binarity of Be stars (see references above), and observation- 
ally this was a completely unconstrained problem. If true, the large 
fraction of (sub-)arcsecond Be star binaries, would have important 
implications for our understanding of Be stars. Obviously, if wide 
binarity could be linked to the Be star phenomenon it would require 
new physics or new models as the effect from distant companions 
will be subtle. Alternatively, low number statistics could have mim- 
icked a high binary fraction, whereas, in reality, this may not be the 
case. A follow-on study with more objects matched by a sample of 
normal B stars should allow us to decide on the issue. Whereas the 
technique of spectro-astrometry is very powerful in detecting bina- 
ries amongst emission-line stars, it is not as sensitive whe n consid- 
ering normal, absorption line, B-type stars (Baines et al. 20o3). In 
order to do a comparative study, we need to employ a technique 
that is both powerful and applicable to both types of star. 

Here we present observations for a simple experiment: are Be 
stars more likely than normal B stars to be in a binary system? 
To this end, we obtained Adaptive Optics (AO) data of a matched 
sample of 40 B and Be stars using an 8 meter class telescope. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target selection 

In order to arrive at an unbi ased, bright sample of Be sta rs, we used 
the Bright Star Catalogue faoffleit & Warren Jn[l98l BSC) as a 
starting point. We selected on spectral type and luminosity class 
(BlVe and BVe), and on Declinations that can be covered by the 
ESO-VLT. There are 71 BVe and BlVe stars in the BSC and its 
Supplement located below +10° Declination. As comparison, we 
selected all BV and BIV stars with Declination smaller than +10°. 
This returned 640 objects, consistent with the notion mentioned 
above that about 15% of all B stars are of Be spectral type. A ran- 
dom sub-sample of B stars was selected from this master sample 
to match the number of Be stars to be observed. As Be stars pre- 
dominately occur for earlier B spectral types, a slight preference 
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Figure 2. The number of 2MASS sources per square arcminute around the 
target objects. A density of 20 arcmin~ 2 corresponds to an average of one 
background object per field of view of 13. 6x 13.6 arcsec 2 . All targets in 
fields with source densities exceeding this value have indeed more than one 
object in the FoV and were initially flagged as having a binary companion, 
but are not included in the final sample. 



was given to the earlier type B stars in this selection. The typical 
K band magnitudes of the Be stars are around 5-7. The intrinsic K 
magnitudes for M0 dwarfs are (only) 6-9 magnitudes fainter than 
those of B9V and B0V stars respectively. Thus, for a limiting K- 
band magnitude of 12.5 almost the entire main sequence down to 
the M dwarfs can be probed. A total of 40 B and 39 Be stars was 
observed using adaptive optics (AO) on Yepun, one of UTs of the 
Very Large Telescope (VLT) at the European Southern Observatory 
(ESO) in Chile. No selection on binarity or brightness was applied, 
instead, the targets were chosen for their observability in terms of 
their Right Ascension, their location on the sky. The targets are pre- 
sented in Table Q] 



2.2 Observations 

The observations were carried out in service mode during 9 nights 
between 24 March 2005 and 5 June 2005 . The observations 
were conducted wi th the AO system NACO dlxnzen et alj|2003 : 
iRousset et alj2003h . In order to achieve the best resolution images, 
data were taken in the K band, where the diffraction limit is 0.057 
arcsec or 57 milli-arcsec (mas). NACO has a 1024x 1024 pixel ar- 
ray, with a pixel size of 27 /im. The finest plate scale of 13.27 
mas/pix was used, giving a field of view of 13.6 arcsec x 13.6 
arcsec. 

As the targets are very bright, the NB2.17 narrow band filter 
was used. This filter is centered on Br7 with a central wavelength 
of 2.166 fim and a full-width at half maximum (FWHM) of 0.023 
/im. The contribution of Br7 emission to the total flux from the Be 
stars is negligible - a typical emission line wit h an Equivalent Width 
of 2-3 A (see e.g. Pers son & McGregoJ 19851) would contribute less 
than 1 % to the total flux received by the filter. For each target star, 
5 frames were taken in the auto-jitter mode with a box size of 5 x 
5 arcsec. The integration time was dependent on the brightness of 
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Figure 3. The binary separations (left) and their cumulative distributions (right) of the B and Be binaries respectively. The Be star data are shown with red 
solid lines, those of the B stars as black dashed lines. Following KS statistics, it can be concluded that the samples are drawn from the same parent distribution 
(see text). 



the target, and was in most cases 12 seconds per jitter position. For 
the very brightest sources, a Neutral Density filter was employed. 
The seeing was always better than 1.4 arcsec, and airmasses were 
below 1.2, in order to achieve a good AO correction. The targets 
were used as AO reference stars. 

The data were reduced with the ESO pipeline reduction soft- 
ware. Dark frames were produced on the same night as each ob- 
serving run, with the same integration time as the science frames 
produced. Dark subtraction and flat field division were performed 
on individual frames by the pipeline, before the 5 calibrated frames 
were combined. This was performed by automated detection of the 
target, determining an offset then registering and stacking of the 5 
calibrated frames. All the raw data were checked by eye to inves- 
tigate whether the automated detection did not produce artificial 
binaries, and in the few ambiguous cases, the shift-and-add proce- 
dure was done by hand. 

Typical FWHM of 0.09 arcsec were measured on the final im- 
ages, and the point spread function's (PSF) widths ranged between 
nightly averages of 0.07 arcsec and 0.10 arcsec. The typical Strehl 
ratio achieved was 30%, with a fair fraction of targets having Strehl 
ratios exceeding 45%. 



2.3 Binary Detection 

The raw and reduced images were visually inspected to identify bi- 
nary companions. The binary separation and position angle were 
determined by fitting a two-dimensional Gaussian function to both 
spatial profiles. The maximum separation found was 8 arcsec, 
which is about the maximum possible value in our square field- 
of-view with the target source positioned in the centre of the frame 
and a 5 arcsec jitter offset. The minimum separation was 0.1 arc- 
sec, at the resolution limit of our data. The difference in A"-band 
magnitude (or rather, difference in magnitude of the light through 



the NB filter, Am, which is very close to AK), was measured for 
each component using either aperture photometry, or by estimat- 
ing the total counts of each star, based on a two-dimensional Gauss 
fit of the spatial profiles. In the case of more than one companion, 
the magnitude difference was found for the closest of the compan- 
ions. In two cases it proved difficult to obtain a reliable estimate 
of the secondary's magnitude and the entry is left blank in Table|2] 
The largest, reliable, magnitude difference was 7.3 magnitude. As 
a check on the quality of our data, the parameters of known bi- 
naries were compared to our results. For example, HR 3745 was 
id entified as having a se paration of 0.55 arcsec and a PA of 275° 
bv lHartkopf et al. ( 1996). Here we find a separation of 0.55 arcsec. 
and a PA of 281°. The difference in position angle is not neces- 
sarily an indication of the errors involved as it can well be due to 
orbital motion. 

We find 14 visual binaries among the 40 observed B stars and 
13 visual binaries among the 39 Be stars observed. Taken at face 
value, the binary fractions are therefore 35 ± 8 % (70% confidence 
interval) for the B stars and 33 ± 8 % for Be stars. 

The numbers for B and Be stars are very close and comfort- 
ably within the uncertainties. However, it is not obvious that these 
concern physical binaries or visual binaries whose secondaries have 
no physical association with the target stars. The best way to test 
this is to obtain second epoch data to investigate whether the sec- 
ondary stars share a common proper motion or not. In the absence 
of such data, we consider the following as a proxy for the possibil- 
ity of finding a chance binary companion. 

The chances of finding one star in the field of view are large 
when the stellar density is at, or exceeds, a value of 1 star per field 
of view of 13.6 x 13.6 arcsec 2 corresponding to 19.5 stars per 
arcmin 2 . A measure of the local surface densit y of stars was ob- 
tained from the 2MASS Point Source Catalogue dCutri et aljioolh 
by tallying up all objects within a circle with a radius of 1 arcmin 
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Figure 4. The masses of Main Sequence stars (black dots) and giant stars 
(red dots) plotted against their intrinsic K band magnitudes. There is a one- 
to-one relation between these parameters for Main Sequence stars, allowing 
the mass ratio of a binary to be derived based on the K band magnitude 
difference only. As we know the spectral type of the primary, we can in 
most cases rule out the secondary to be a giant or supergiant, as these would 
be brighter than the primary, which is not observed. 



from the target objects. The frequency of the stellar densities to- 
wards the sample is shown in Fig. [2] Most objects lie in regions 
with stellar densities less than 10 per arcmin 2 , less than half a star 
per field of view, and the presence of a chance companion is not 
expected. However, all six objects which lie in densely populated 
areas on the sky (densities larger than 20 arcmin -2 ), are found to 
have a companion. It is much more likely that these objects are 
normal field stars that happen to be in the field of view rather than 
physical binary stars, and we will proceed without those objects. 



3 RESULTS 

3.1 Binary fraction 

After selecting on the objects with a local field star surface den- 
sity less than 20 arcmin" 2 , we remain with 11 Be binaries out of 
37 targets (30 ± 8 %) and 10 B binaries out of 36 targets (29 ± 
8%). The B and Be binary fractions are remarkably similar. In the 
remainder of this study, we will continue the analysis with these 21 
binary systems. For completeness we note that if the density cri- 
terion were lowered to 12.5 or 15 to exclude the next two objects, 
the statistics remain the same, but now both samples would be re- 
duced by exactly one object. HR 5661 (a Be star with a local source 
density of 16 arcmin -2 ) and HR 6316 (a B star with a local source 
density of 18 arcmin -2 ) remain in the sample. 

In the following, we investigate and compare basic properties 
of the binary systems such as separations and flux differences - 
which directly translate into mass ratios - between the primaries 
and secondaries. The binary systems and their parameters are listed 
in Table [2] 



3.2 Binary separation 

Here we compare the binary separations for all B and Be stars. A 
histogram of the separations is shown in Fig. [3] The Be-star sep- 
arations are slightly more centrally peaked than the B star separa- 
tions. Although it is not a very significant result, we would wish to 
draw attention to the fact that this may be expected. The Be stars 
are rarer than B stars and they are thus, on average, further away. 
This is corroborated by the parallaxes of th e sample objects . Using 
th e Hipparcos parallaxe s taken either from lPerrvman et sil d 19971) 
or Ivan Leeu wen (2007), we find that the average parallax of the 
B-binaries is 4.9 milli-arcsec (with a generous spread of 2.7 mas) 
while the average parallax of the Be stars is slightly lower at 4.4 
mas (with a a of 2.8 mas). Should the intrinsic separation distribu- 
tion be similar, then it is to be expected that the Be stars' separations 
(in arcsec) are smaller, or somewhat more centrally peaked than for 
the B stars. We note however, that this is not a large effect. We also 
performed a 2-sample Kolmogorov-Smirnov (KS) test. The cumu- 
lative distribution is shown in the right hand panel of Fig. [3] The 
KS statistic indicates that the hypothesis that the Be and B binary 
separations are drawn from different separation distributions can 
be rejected at the 74% level. This is not overly conclusive, but the 
sample is large enough to have identified large differences in the 
distribution of separations, if there would have been any. 

As conclusion of this section, we find that the Be binaries are 
similar to the B binaries in terms of separation. 



3.3 Mass ratios 

Here we consider the masses of the components in the binary sys- 
tems. In principle, we can infer the masses of the primaries from 
their listed spectral types using published calibrations, and derive 
spectral types of the secondaries using the flux differences and 
by comparing with the intrinsic magnitudes listed for the differ- 
ent spectral types, we can then determine the secondaries' masses. 
Although this is a standard method, it is fairly elaborate and intro- 
duces various errors. These comprise the error in spectral typing, 
the mass and magnitudes range f or a given spectral typ e (cf. e.g. 
the width of the Main Sequence, IPerrvman et al.1 dl997l) . see also 
Oudmaijer et al. (1999) for an in-depth study of KOV stars) and af- 
fect the final result. 

These issues are mostly avoided in the present study by sim- 
ply using the K band magnitude difference. This provides an ex- 
cellent opportunity to derive the mass ratio with no dependence 
at all on the spectral type of either component. It turns out that 
there is a strong relation between the intrinsic K band magni- 
tude and the mass of a main sequence st ar. This is shown in 
Fig- El where the mass of a star (taken from lStraizvs & Kurilienel 
1981) is plotted against the intrinsic K band magnitude which was 
computed using the i ntrins ic V band magnitudes also taken from 
Straizys & Kurilienel dl98ll) and the intrinsic V — K colours due to 



Koornneef Jl983l) . The relation for Main Sequence stars is close to 



a straight line which can be represented as 



Iog(M) 



-0.18^ + 0.64 



with M in solar units, and a formal uncertainty on the slope of 
0.005. This in turn translates into a mass ratio q - without the need 
of knowing the masses separately: 



10" 
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Figure 5. As Fig. [3] but now for the mass ratios of the binary systems. The mass ratios of the B and Be star binaries are no t markedly different and acco rding 
to KS statistics drawn from the same parent population. The dotted line is not a fit to the data, but the distribution found by Shatsky & Tokovinin 1 2002) for a 
larger sample of B stars. 



with AK the magnitude difference between both objects in our 
data. 

By virtue of the rather linear relationship between K band 
flux and stellar mass, the mass ratio can be determined indepen- 
dent of prior knowledge of the stellar spectral types and luminosi- 
ties while it is independent of reddening too. The only assump- 
tion we tacitly make is that both the primary and secondary stars 
under consideration are on the Main Sequence. This may seem 
obvious, but we note that many known binaries with a hot O, B 
typ e dwarf star have a cool, G, K or M-type giant secondary (see 
e.g. iGinestet & Carquiliatll2002h . In our sample, HR 2577, an un- 
resolved spectroscopic binary with a K2II companion listed in Ta- 
ble [T] is a prime example. In most such binary stars, the B-type 
component dominates at optical wavelengths, while the secondary 
dominates at longer wavelengths, simply because cooler stars are 
redder. In contrast, the current study is conducted in the infrared al- 
ready and the brightest binary component is the previously known 
B-star whose catalogued A* band magnitide is consistent with the 
optical brightness, while the secondary is fainter. 

For reference, we have also plotted the mass of a selected 
number of giant stars against their intrinsic K band magnitudes in 
Fig. [4] The cool G, K and M-giants are brighter than most B stars 
up to a spectral type of B3V. This means that for all stars with spec- 
tral types of B3 or later in our sample, we can be fairly certain that 
the companion is not a giant as the primary star, the B star, is the 
brightest. Indeed, only for B0V-B3V stars there would be a chance 
that the (fainter) companion is a giant rather than a main sequence 
dwarf star. In such a case, the magnitude difference will be at most 
2 magnitudes. Inspection of Table [2] however reveals that almost 
all binaries with a B0V-B3V primary have a magnitude difference 
larger than 2. We conclude that it is very likely that the secondary 
stars are main sequence stars and that we are therefore justified in 



applying our Admass relation to determine the mass ratio of the 
systems. 

One correction needs to be made for the Be stars, the free-free 
emission from the ionized disk increases with wa velength, and can 
give a slight excess emission at the A"-band (e.g. iDoughertv et al.l 
1 199 lh . To arrive at the intrinsic magnitude of the star itself, we 
need to subtract the disk contribution from the K band flux. This 
is done in the standard manner by dereddening the V band magni- 
tude (using intrin sic colours for the respective spectral types from 
Schmidt-Kaler in lAller et alj fl982l Chap. 4)) and computing the 
predicted K band magnitude using the V — K colours as outlined 
above. The difference between the predicted and observed A"-band 
magnitude is the excess. In two cases, HR 6712 and HR 5193, the 
excesses were negative by up to half a magnitude and thus unphys- 
ical. These stars are variable and using the optical colours from 
SIMBAD rather than from the BSC solved this situation. The "ex- 
cesses" computed for the normal B stars averaged 0.01 magnitude 
with a scatter of 0.1 magnitude, which gives an indication of the 
uncertainties involved in this exercise. As all magnitude differences 
between the companions are much larger than 0.3 magnitudes, this 
error does not affect the resulting mass ratio. In fact, only two Be 
stars have computed excesses which are much larger than 0.3 mag- 
nitude, these are HR 2142 (0.6 mag) and HR 2356 (0.8 mag). After 
correcting for the excess emission, we can then compute the pho- 
tospheric magnitude difference and the mass ratio q. 

In Fig.|5]the results for the mass ratios are presented in a sim- 
ilar fashion as for the separations in the previous figure. Note that 
the two objects for whose close companions we were unable to re- 
trieve accurate photometry are not plotted. The left hand plot shows 
the distribution of mass ratios, which range from 0.07 to 0.9. Based 
on the relatively low numbers, it is hard to describe this distribu- 
tion, but it would seem that the distributions are similar and rather 
flat for the B and Be stars. It does not rise steeply to low masses as 
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Table 2. The final binary detections, with objects in dense fields excluded 
(see text). Listed are the separations and position angles (measured East of 
North), as well as the magnitude differences. Typical errors are of order 0.05 
arcsec, 1° and 0.1 magnitude respectively. K band magnitudes are taken 
from 2MASS. The mass ratios q are derived in the text. 

would be the case if the secondary mass es were randomly sample d 
from the Initial Mass Function (see e.g. IWheelwright et al.ll2010h . 
The cumulative plot in the right hand panel confirms that there is 
not a significant difference between the two distributions. The KS 
test returns a probability that the samples are drawn from differ- 
ent populations can be rejected at 50%. As with the separations, if 
the mass ratios would have been very different, the current sam- 
ple is large enough to highlight this. It therefore seems likely that 
the mass ratios of both the B and Be type binaries too are sam- 
pled from the same parent distribution. For reference the function 
f(q) oc q~ 5 is overpl otted as a dashed line. This is not a fit to our 
data, instead it is what IShatskv & Tokovinirj ( 2002) found to best 
represent the mass ratio they determined for a large sample of B- 
type binaries. Our data are consistent with their data, which sample 
a similar separation and magnitude range. 

In summary of this exercise, there is no difference between B 
and Be binaries as far as their companion masses is concerned. 



4 DISCUSSION AND FINAL REMARKS 

We have observed a total of 79 B and Be stars at high spatial res- 
olution to study their binary properties. This is the first system- 
atic comparative and deepest imaging study into the binarity of Be 
stars published. The observations probe scales down to less than 
0.1 arcsecond and allow companions up to 10 magnitudes fainter 
to be detected. In the process, we derived a simple relation between 
the A"-band magnitude difference and the mass-ratio of the binary 
components under the assumption that both components are on the 
Main Sequence. The latter is justified in our sample. After exclud- 
ing a few visual binaries that are located in regions of high stellar 
density, we find that : 



(i) There are 10 detected binaries out of 36 B stars and 11 de- 
tected binaries out of 37 Be stars (corresponding to binary fractions 
of 29 ± 8 % and 30 ± 8% respectively). The binary fractions are 
nearly identical. 

(ii) The separations of the B and Be binary systems follow the 
same cumulative distributions. 

(iii) The mass ratios for the 9 B and 10 Be binaries for which we 
have differential photometry have similar distributions. Like for the 
separations, the samples appear to be drawn from the same parent 
distribution. 

This leads to our main conclusion : 

• The binary fraction, and the binaries' properties, of Be stars 
in the separation range 20 - 1000 AU are very similar to those of 
a matched sample of B stars. This is a statistically robust result 
and as both samples were observed with the same observational 
equipment, and therefore were subject to the same observational 
biases, we can rule out binary companions at separations of 20 AU 
or more as being responsible for the Be phenomenon. 

Returning to the original motivation of the survey, the previ- 
ously found large fraction of wide Be binaries has been reduced 
by a factor of 2 due to improved number statistics. More impor- 
tantly, the comparison with normal B stars now allows us to firmly 
conclude that there is no need, based on this study, to explore new 
models involving wide binaries for the Be phenomenon. 

Closer binaries than 20 AU can only be found with methods 
that have even higher spatial resolution than here (diffraction lim- 
ited imaging at an 8 meter telescope in the A* band), or with radial 
velocity studies. A systematic imaging study at even higher res- 
olution (indeed at any resolution) such as this one has not been 
perfo rmed, but the systematic spectroscopic study by dAbt & Levy! 
Il978l) does not reveal a large percentage of close binaries. Indeed, 
on the contrary, they find a lower close binary fraction among Be 
stars than B stars. Therefore, binarity can not be responsible for all 
Be stars. 

Finally, a word on the similar mass ratios and binary separa- 
tions. The mass ratio and separation of a binary system are the re- 
sult of the initial conditions governing the formation of both com- 
ponents and their subsequent evolution. The stars in our sample 
are field stars, which means that if they were formed in clusters, 
as most massive stars do, these clusters have dispersed long since. 
Given that the separations of the B and Be binaries under consid- 
eration are fairly wide, the objects were most likely formed in low 
density clusters or even in isolation. This is because large separa- 
tion binaries will be destroyed early in the evolution of dense clus- 
ters, while close, more strongly bound binaries are most likely to 
survive throughou t the entire evolution jKouwenhoven et al]|2009l : 
iParker et al. 2009). It is beyond the scope of this paper to derive the 
precise cluster conditions, but it would appear that both the evolu- 
tionary history and the birthplaces of the B and Be stars are very 
much the same. It can be concluded from our data, therefore, that 
the Be phenomenon does not result from a substantially different 
star formation mechanism or conditions compared to those for nor- 
mal B stars. 
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